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INTRODUCTION 


Supersonic  molecular  jet  expansions  are  known  to  produce  a  rich 

array  of  van  der  Waals  (vdW)  clusters  in  an  environment  accessible  to 

various  spectroscopic  techniques.^  Photoionization  time  of  flight 

mass  spectroscopy  (TOFMS)  is  especially  useful  for  the  deconvolution 

2 

of  the  absorption  spectra  of  various  clusters  in  the  beam  .  The 

study  of  large  clusters  can  be  carried  out  with  this  relatively  easy 

mass  selective  technique  and  thereby  the  modeling  of  condensed  phases 

and  aggregated  systems  can  be  explored.  Furthermore,  the  evolution  of 

the  mass  detected  absorption  spectrum  of  a  central  chromophore  as 

larger  clusters  are  observed  provides  some  insight  into  nucleation 

processes  and  intermolecular  interactions  in  condensed  phases. 

One  readily  attainable  spectroscopic  parameter  which  can  be  of 

help  in  characterizing  a  cluster  is  the  spectral  shift  for  a  given 

transition;  that  is,  the  difference  in  energy  between  a  cluster 

vibronic  feature  and  the  comparable  feature  in  the  isolated  solute 

molecule.  It  is  found  that  in  the  case  of  benzene  solvated  with 

methane,  ethane,  or  propane,  the  spectral  shift  depends  essentially  on 

the  cluster  geometry.  The  general  geometries  of  most  of  the  small 

clusters  can  be  deduced  from  the  qualitative  analysis  of  the  patterns 

and  trends  of  the  spectral  shifts,  from  vdW  vibronic  features,  and  by 

•5 

analogy  with  other  systems  for  which  the  geometry  is  known.  For  many 
clusters,  additional  information  is  necessary  to  deduce  cluster 
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geometry.  One  approach  is  to  use  high  resolution  spectroscopy  to 
obtain  the  rotational  spectra  of  these  clusters. Aside  from 
being  technically  difficult,  this  approach  usually  will  not  yield 
detailed  enough  information  to  reveal  the  exact  structures  of  larger 
clusters.  A  simpler  and  potentially  more  useful  approach  to  finding 
the  geometries  of  larger  clusters  is  to  use  computer  modeling  of  the 
clusters  to  calculate  minimum  energy  configurations.  If  the  low 
energy  configurations  of  a  given  cluster  are  consistent  with 
spectroscopic  data  for  the  system,  the  spectroscopic  features  of  a 
particular  geometry  can  be  identified.  The  assignments  help  to 
elucidate  the  complicated  relationship  between  geometry  and  spectral 
shift. 

The  assignment  of  spectral  features  to  specific  geometries  also 
provides  insight  into  the  nucleation  process  in  the  jet.  Solvent 
molecules  may  add  to  benzene  one  molecule  per  collision,  a  process  we 
will  call  homogeneous  nucleation,  or  a  benzene  molecule 
may  bind  to  a  solvent  cluster  to  add  more  than  one  molecule  per 
collision.  This  latter  inhomogeneous  nucleation  typically 
results  in  very  anisotropic  clusters  -  clusters  in  which  most  of  the 
solvent  molecuxes  are  on  one  side  of  the  ring.  Homogeneous  nucleation 
can  produce  relatively  isotropic  clusters  which  m  ly  be  thought  of  as 
model  liquid  solution  or  matrix  systems.^  However,  inhomogeneous 
nucleation  seems  to  predominate  in  the  formation  of  the  larger 
clusters  in  this  study. 
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This  report  will  present  the  spectroscopic  results  of  the  study 
of  benzene  clustered  with  methane,  ethane  and  propane  (Ben(CH  )  , 

^  X 

Ben(C2  Ben(Cj  together  with  the  results  of  computer 

modeling  of  the  clusters.  Attention  will  be  confined  in  this  paper  to 
small  clusters  (x<3)*  In  a  companion  paper  (hereafter  referred 
to  as  TI)  the  results  for  toluene  clustered  with  methane,  ethane,  and 
propane  will  be  presented.  The  differences  and  similarities  between 
benzene  clusters  and  toluene  clusters  will  be  discussed  in  TI. 

II  EXPERIMENTAL  PROCEDURES 

The  apparatus  and  techniques  employed  in  this  study  are  largely 

Q 

the  same  as  those  previously  reported.  Fluorescence  excitation 

experiments  were  found  to  be  of  limited  utility  as  the  clusters  of 

interest  are  nearly  always  minority  species  in  the  beam.  In 

addition,  dispersed  emission  from  most  of  the  clusters  is  weak  and 

g 

complicated  as  noted  previously.  Consequently,  only  mass  selective 
absorption  spectra  (2-color  TOFMS)  will  be  discussed  in  this 
paper. 

The  optimum  nozzle  backing  pressure  is  the  same  for  all 
spectra  reported  in  this  work  (P^  '^'100  psi)  and  the  nozzle  is  kept  at 
ca.  300K  throughout  these  experiments.  The  concentration  of  solvent  in 
the  backing  region  for  all  systems  discussed  is  2%.  Dispite  nearly 
identical  backing  conditions  throughout,  intensity  comparisons  between 
species  are  difficult.  The  various  solvents  cluster  differently  and  a 
particular  expansion  condition  can  produce  different  concentrations  of 
multimers  for  different  solvents. 
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Otner  changes  have  been  made  in  the  previously  reported 
experimental  conditions.  A  Comlinear  Corp  CLC  1UO  linear  amplifier  is 
placed  directly  after  the  ion  detector  as  required  to  enhance  the 
weakest  signals  by  a  factor  of  10.  The  viton  0-ring  which  seals  the 
pulsed  nozzle  has  been  replaced  by  a  Teflon  coated  viton  0-ring  to 
present  swelling  of  the  0-ring  due  to  absorption  of  benzene. 

Finally,  a  cost  savings  of  about  $20/10  hours  of  laser 
operation  time  can  be  realized  by  recycling  the  coumarin  bOO  laser 
dye.  Most  of  the  undecomposed  dye  can  be  recoverc;d  from  the  batch  of 
used  dye  (which  lasts  about  20  nours)  by  a  single  recrystallization 
from  methanol  and  water.  31ow  evaporation  of  a  dye-methanol  solution 
produces  good  quality  dye  crystals:  no  reduction  in  performance  is 
noted  for  recycled  dye. 

The  calculation  of  cluster  energy  is  performed  on  an  HP  9845 
desktop  computer.  An  atom  -  atom  (exp-bj  intermolecular  potential  of 
the  form 


E  .  =  -  — expi-C  R  .1 

ij  „  6  ij  ij  ij-* 


IS  employed  in  which  A,H  and  C  are  parameters  and  is  the 
distance  between  the  i^^  and  atoms  of  different  molecules. 

Summing  over  all  pairs  of  atoms  on  different  molecules  gives  the  total 
intermolecular  energy  of  the  cluster. 

Parameters  used  in  these  calculations  are  adapted  from  those 
given  in  the  literature. ^ ^ ^  Parameters  determined  for  methane  are 
used  for  all  of  the  solvents.  The  values  used  are  given  in  Table  I. 
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The  minimization  routine  to  find  the  minimum  energy 
configurations  involves  moving  one  solvent  molecule  at  a  time  a  step 
in  each  of  the  translational  and  each  of  the  three  rotational  degrees 
of  freedom.  The  new  energy  of  the  cluster  is  calculated  and  compared 
to  the  previous  energy.  If  the  new  energy  is  lower,  the  step  is 
repeated;  if  it  is  not  lower,  the  next  degree  of  freedom  is  tried. 
Although  simple,  this  technique  cannot  guarantee  that  the  absolute 
minimum  energy  configuration  has  been  found.  Consequently,  local 
minima  in  the  potential  energy  surface  are  often  found  and  one  must  be 
careful  to  sample  all  of  configuration  space  to  be  certain  that  the 
lowest  energy  configuration  has  indeed  been  located.  Often  intuition 
developed  through  the  study  of  simpler  systems  can  limit  the  number  of 
possible  initial  configurations  that  must  be  sampled. 

The  calculated  binding  energy  is  the  depth  of  the  ground  state 

I 

vdW  potential  well  (D^").  The  measured  binding  energies  are  . 

The  difference  between  these  values  includes  ground  and  excited  state 
energy  differences  (cluster  spectral  shifts)  and  zero  point  energies. 
Considering  the  accuracies  of  the  experiments  and  the  calculations, 
these  differences  are  at  present  neglected  in  the  discussion. 

Ill  RESULTS 

The  allowed  6^  transition  of  benzene  is  intense  in  all  of  the 
0 

clusters  and  is  relatively  free  of  influence  from  other  transitions. 
The  spectra  of  the  various  clusters  are  easily  assigned  in  this 
region.  The  0^  region  was  also  scanned  for  each  of  the  clusters: 
only  clusters  with  symmetry  lower  than  C^  will  exhibit  a  0^  origin 
peak.  The  presence  or  absence  of  the  0^  feature  of  a  cluster  helps 
to  identify  the  cluster  geometry.  While  vibronic  transitions  higher 


than  of  benzene  are  observed,  no  higher  transitions  of  the 
clusters  are  reported  because  of  congestion  and  low  intensity  of  the 
cluster  spectra. 

The  assignments  for  vdW  vibrational  features  given  below  are 
made  to  be  consistent  with  previous  assignments  of  the  aniline- 
methane,^  and  toluene-methane'^  systems  and  with  assignments  given  in 
II.  In  all  of  these  systems  the  frequencies  of  the  symm.  trie  vai«/ 
stretches  scale  appropriately  with  the  masses  and  force  constants 
(binding  energies).  The  vdW  bending  modes  are  typically  lower  in 
energy  and  intensity  than  the  stretches.  Using  tnese  general 
guidelines  the  assignments  given  below  are  reasonable  and  fit  well 
into  the  previous  patterns  of  vdirf  mode  frequencies. 

Benzene-Me thane 

The  b^  transition  of  Ben(CH^)^  is  shifted  -41cm~^  from 
the  benzene  b^  (see  figure  1  and  Table  II).  A  progression  in  the 
VdW  stretch  is  observed  with  a  stretching  frequency  of  5i£  cm  ^  .  In 
addition  a  vdW  bend  with  a  frequency  of  27  cm  ^  is  seen.  The  0^ 
of  Ben(CH^)^  is  not  observed,  indicating  that  the  5-fold  axis  of 
benzene  is  preserved  in  tne  Ben(CH^;^  cluster. 

Ben(Uri^)2  at  b^  reveals  two  configurations  (see  figure  1, 
Table  11) •  Configuration  1  has  a  spectral  shift  of  -b1 cm  S  twice 
that  of  Bt,n(CH^)p  A  symmetric  stretch  with  a  frequency  of  51  cm  ^ 

13  noted.  Configuration  II  of  Ben(CH^)^  evidences  a  -40  cm  ^ 
spectral  shift:  this  feature  is  broader  than  those  of  configuration 
i . 
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In  the  regxon  of  the  Ben(CH^;^  0^  only  one  peak  with  a 
spectral  shift  of  -40  cm  ^  is  observed  (.see  figure  2):  therefore, 
configuration  I  of  Ben(CH^J2  preserves  the  :5-fold  axis  of  benzene 
and  configuration  II  does  not. 

The  geometries  of  Ben(,CH^)^  2  to  a  large  extent,  be 

deduced  from  the  spectra.  In  Ben(.CH^)^  the  methane  must  rest 
symmetrically  above  the  ring  to  preserve  the  '.5-fold  axis  of  benzene. 
Configuration  I  of  BenCCH^j^  must  nave  methane  symmetrically  above 
and  below  the  benzene  ring,  ana  configuration  II  must  have  both 
methanes  on  the  same  side  of  the  ring  to  produce  an  asymmetric 
geometry  with  a  non-additive  shift  and  an  induced  transition. 

The  computer  moaeling  of  these  two  clusters  agrees  perfectly 
with  the  geometries  deduced  from  the  spectra  (.see  figure  5,  Table 
lllj.  The  mj.nimium  energy  configuration  of  Ben(,CH^).|  preserves  the 
j>-fold  axis  of  benzene.  The  predicted  binding  energy  is  599  im  S  a 
value  greater  than  5^2  cm  ^as  required  if  the  is  to  be  observed 
for  BenvCH^j^.  Two  configurations  for  Ben^Cri^j^  arise  from  the 
computer  modeling.  One  preserves  the  5-fold  axis  of  benzene  and 
should  have  a  spectral  shift  twice  that  of  Ben(,CH^)^,  and  the  other 
configuration  is  asymmetric  with  both  methanes  on  the  same  side  of  the 
ring. 

The  spectrum  of  Ben(CH^).^  evidences  two  configurations  also. 
Configuration  I  has  a  shift  similar  to  configuration  I  of  Ben(,CH^J2 
(.see  figure  1j.  This  geometry  probably  consists  of  Ben(.CH^;^ 
configuration  II  with  a  methane  bound  symmetrically  below  the  ring  for 


a  total  shift  of  about  -80  cm  ^ .  The  other  configuration  of 
Beni,CH^)^  probably  features  all  three  methanes  on  the  seme  side  of 
the  ring  to  produce  a  small  (-40  cm  spectral  shift.  The  peaks  in 
the  Ben(CH^)^  spectrum  are  more  diffuse  than  those  of  Ben 
(^^4^1  2  indicating  that  a  number  of  slightly  different 
configurations  probably  exist.  The  Ben(CH4)^  fa^  spectrum  is 
sufficiently  weak  that  the  0^  is  not  to  be  expected.  Calculations 
were  not  done  for  this  cluster. 

Benzene-Ethane 

BenCC^H^)^  shows  two  different  configurations  in  the  6^ 
region  (see  figure  4,  Table  IV).  Configuration  I  has  a  spectral  shift 
of  -'ol  cm  ^  and  a  long,  regular  progression  in  .he  vdW  stretch  with 
a  frequency  of  17  cm  Configuration  II  has  a  shift  of  -tjl  cm~^ 
and  a  stretching  frequency  of  16  cm”\  The  spectrum  of 
Ben(C2H^)^  at  the  0^  (see  figure  b.  Table  V)  exhibits  only  the 
origin  associated  with  configuration  II  and  two  associated  vdW  bends. 
The  relative  intensity  of  the  bends  is  much  larger  at  the  origin  tnan 
at  6^;  the  vdW  stretch  is  not  evident  at  the  origin,  perhaps  due  to 
the  similarity  of  the  ground  and  excited  state  potentials. 

The  Ben(C2H^)2  (see  figure  4,  Table  11)  configuration  I  has 
a  spectral  shift  of  cm  ' ,  at  6^  and  is  not  seen  at  the  0^. 
Configuration  II  of  this  species  reveals  a  spectral  shift  of  -4b 
cm  \  a  stretching  fiequency  of  1b  cm  \  and  one  bending  mode  in 
the  by  region.  At  the  origin,  only  the  of  configuration  II  is 
found  with  Its  associaied  vdW  bends  (see  figure  b,  Table  V);  again, 
no  stretches  are  observed. 


The  calculations  for  Ben(C„H,  ),  agree  quite  nicely  with  the 

<1  D  1 

spectra  (see  figure  6,  Table  VI Two  configurations  of  minimum 
energy  are  found:  one  which  preserves  the  3-fold  axis  of  benzene  and 
one  which  does  not.  The  asymmeiric  configuration  is  predicted  to  have 
a  1 20  cm  ^  larger  binding  energy. 

The  computer  modeling  of  Ben(C2H^)2  reveals  several  stable 
geometries  two  of  which  are  presented  in  figure  7  and  Table  VII.  All 
geometries  involving  an  ethane  on  either  side  of  the  ring  must  be 
dismissed  since  the  spectrum  of  Ben(C2H^^2  region  of  6^ 

shows  no  configurations  with  an  additive  spectral  shift  with  respect 
to  Ben(C2H^)^.  Only  one  geometry  (labeled  1  in  figure  7)  is  found 
in  the  computer  modeling  which  preserves  the  '3-fold  axis  of  benzene. 
Several  other  geometries  of  Ben(C2H^J2*  which  should  have  spectral 
shifts  similar  to  that  found  for  configuration  II  (see  figure  4),  are 
predicted  in  the  computer  modeling  of  Ben(C2H^)2'  The  most  stable 
of  these  is  shown  in  figure  7- 
Benzene -Pro pane 

Oniy  one  configuration  of  Ben(C^Hyj^  is  clearly  evident  from 
the  spectrum  (figure  b,  Table  Vlli^.  The  peak  of  this  cluster  is 
shifted  -72  cm  ^  from  the  monomer  b^.  An  intense  progression  in 
the  stretch  is  evident  witn  a  frequency  of  24  cm  ^ .  The  unassigned 
peak  -j>7  cm  '  from  the  monomer  b^  may  be  due  to  a  different 
configuration,  or  it  may  be  due  to  some  high  frequency  vdW  motion.  Wo 
0^  peak  13  obse.ved  for  this  cluster  indicating  tnat  the  cluster 
geometry  preserves  the  3-fold  axis  of  benzene. 


The  vdW  stretch  of  BemC^Hy)^  seems  to  have  a  positive 
anharmonicity  for  the  first  two  quanta  of  the  stretch.  Underlying 
peaks  due  to  slightly  different  configurations  could  be  modifying  the 
iineshapes  of  some  of  the  peaks  to  produce  this  effect.  It  is 
difficult  to  envision  how  the  symmetric  stretch  of  the  cluster  could 
have  a  positive  anharmonicity;  however,  this  possibiiiiy  cannot  be 
dismissed.  A  similar  but  much  leas  pronounced  effect  is  observed 
for  the  stretch  of  Benl,C^H^)^  . 

There  is  no  peak  in  the  BenCC^Hy)^  spectrum  with  a 
shift  twice  that  of  Ben(.C.^Hy)^ .  One  configuration  of  Ben(C.^Hy)2 
seen  in  the  b^  region  has  a  spectral  shift  of  -8^  cm  ^ ,  and  a 
stretching  frequency  of  25  cm  \  Other  geometries  with  smaller 
spectral  shifts  are  observed.  No  peaks  are  observed  at  the  cluster 
Oq  position. 

The  computer  modeling  of  Ben(C^Hy)^  shows  several  possible 
configurations.  The  lowest  energy  configuration  possesses  the  highest 
symmetry  and  is  pictured  in  figure  9  (see  Table  IX).  The  three 
hydrogen  atoms  closest  to  the  ring  are  positioned  almost  symmetrically 
above  the  ring.  The  0^  transition  of  Ben^C^Hg)^  is  not  observed 
presumably  due  to  the  approximate  three-fold  symmetry  of  this  cluster. 
The  spectrum  of  Ben(C^Hg;2  does  not  show  an  origin  either.  One 
example  of  a  Ben(C^Hg)2  which  preserves  this  local  symmetry 

IS  pictured  in  figure  9*  Other  similar  configurations  are  possible. 

The  spectrum  of  Ben(C^HgJ-^  in  the  region  indicates  that 
at  least  two  configurations  of  this  cluster  exist  (see  figure  10, 

Table  X).  Due  to  the  small  spectral  shifts,  the  clusters  must  involve 
all  5  solvent  molecules  situated  above  the  ring.  The  cluster  0^ 
transitions  with  associated  vdW  bends  are  seen  for  both  of  these 
configurations.  No  calculations  were  performed  on  these  clusters. 

-  10  - 
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IV  DISCUSSION 


The  ability  to  assign  geometries  to  clusters  giving  rise  to 
certain  peaks  in  the  spectra  leads  to  a  greater  understanding  of  both 
the  energetics  of  the  system  and  the  process  of  cluster  formation.  The 
relationships  between  geometry  and  spectral  shift  will  be  explored 
first,  followed  by  a  discussion  of  the  nucleation  processes  in  the 
jet. 

Spectral  Shift 

The  spectral  shift  for  a  cluster  is  the  difference  between  the 
cluster  binding  energy  in  the  ground  and  excited  states.  The  shift 
depends  on  the  polarizability  of  the  solvent  along  the  intermolecular 
bond  and  on  the  change  in  polarizability  of  the  solute  upon  excitation 
in  the  region  (direction)  of  this  bond.  The  solute  polarizability  may 
change  quite  differently  for  different  binding  sites  on  benzene  to 
produce  different  spectral  shifts  for  different  geometries. 

The  spectra  of  Ben(CH^)^  2  (see  figure  1 )  demonstrate  that 
binding  of  a  methane  symmetrically  above  or  below  the  aromatic  ring 
contributes  about  -40  cm  ^  to  the  spectral  shift  of  the  complex. 
However,  the  asymmetric  binding  of  two  methanes  above  the  ring 
(yielding  an  observed  Oq  transition)  produces  a  shift  of  only  -40 
cm  The  computer  modeling  of  the  asymmetric  Ben(CH^)2  cluster 
shows  that  one  of  the  methanes  is  largely  above  the  ring  (although 
displaced  from  the  center  by  .695  ?)  while  the  other  methane  is 
displaced  by  2.778  %  from  the  ring  center.  This  latter  CH^  is  above 
the  midpoint  of  the  C-H  bond  (see  figure  '5  and  Table  i).  If  this 
geometry  is  correct,  almost  all  of  the  spectral  shift  of  this  cluster 
must  be  due  to  the  interaction  of  benzene  with  the  methane  above  the 


ring.  The  other  methane  must  contribute  virtually  nothing  to  the 
spectral  shift  even  though  it  is  bound  nearly  as  tightly  to  benzene. 
In  fact,  the  asymmetric  (II)  BenCCH^)^  cluster  has  a  virtually 
identical  binding  energy,  but  a  much  smaller  spectral  shift,  than  the 
symmetric  BenlCH^)^  cluster. 

These  observations  lead  to  the  conclusion  that  interaction  of 
the  solvent  with  the  aromatic  tr-doud  of  benzene  is  responsible  for 
most  of  the  spectral  shift  even  though  the  binding  energy  may  be 
nearly  as  large  for  other  binding  sites.  Such  a  generalization  helps 
to  explain  the  spectral  shifts  of  other  clusters  as  well 

The  spectral  shifts  of  the  two  distinct  (symmetric  and 
asymmetric)  geometries  of  Ben(C2H^J^  are  deduced  and  assigned  in 
the  Results  section.  The  higher  symmetry  cluster  (figure  b),  which 
preserves  the  three-fold  axis  of  benzene,  has  a  larger  spectral  shift 
than  the  lower  symmetry  cluster  for  which  the  origin  is  observed. 

As  in  the  benzene-methane  system,  the  calculated  binding  energies  for 
the  two  configurations  are  virtually  identical,  yet  the  asymmetric 
configuration  has  a  much  smaller  spectral  shift. 

Apparently,  both  orientation  and  position  of  a  solvent  molecule 
play  an  important  role  in  determining  the  spectral  shift  and  binding 
energy  of  a  cluster  of  given  composition.  Ethane  has  a  larger 

1 2 

polarizability  parallel  to  its  long  axis  than  perpendicular  to  it. 
Therefore  the  orientation  of  the  long  axis  of  ethane  with  respect  to 
the  plane  of  the  benzene  ring,  and  the  position  of  the  center  of  the 
ethane  molecule  with  respect  to  the  center  of  the  benzene  ring  have 
considerable  effect  on  both  spectral  shifts  and  binding  energies.  The 


larger  spectral  shift  is  associated  with  the  solute-solvent  relative 
orientation  for  which  the  direction  of  large  solvent  polarizability  is 
perpendicular  to  the  plane  of  the  ring,  directly  over  the  center  of 
the  molecular  system.  The  asymmetric  configuration  (ll)  of 
Ben(C2H^j^  has  a  larger  calculated  binding  energy  (-77b  vs.  -6154 
cm~^ )  and  a  smaller  shift  (-b1  vs.  -1)7  cm  S  most  likely  arising  from 
a  reduced  iHsyatem  interaction  but  a  larger  overall  benzene-ethane 
interaction:  in  this  configuration  the  ethane  long  axis  is  nearly 
parallel  to  the  plane  of  the  benzene  ring.  Ethane  in  these  two 
different  orientations  interacts  with  benzene  as  though  it  were  two 
different  solvents  with  different  sizes  and  polarizabilities. 

The  spectral  shifts  for  the  two  configurations  of  Ben(C2H^)2 

are  neariy  identical  to  the  shifts  for  the  Ben(C2H^)^ 

configurations.  Based  on  the  calculations,  the  geometries  for 

Ben(C„H,  involve  Ben(C„HA).  structures  with  a  second  CoH, 

C  Q  d  d  o  \  do 

coordinated  on  the  same  side  of  the  ring  as  the  first  C2H^  group: 
the  amount  of  solvent  Tr-cioud  overlap  is  about  the  same  in 

Ben(C2H^)2  Ben(C2H^j^. 

The  calculated  geometry  of  Ben(C.^Hy)^  which  agrees  with  the 
spectroscopic  data  (i.e.  no  0^  transition,  large  spectral  shift) 
features  a  large  interaction  between  the  propane  and  the  ir-system  of 
benzene.  The  calculated  configuration  is  depicted  in  figure  y:  this 
arrangement  of  solute  and  solvent  has  approximate  high  symmetry  with 
all  three  carbon  atoms  lying  near  th  aromatic  n-system.  The  binding 
energies  and  shifts  for  the  Ben(C^Hy)2  clusters  are  also  consistent 
with  the  above  conclusions. 


The  Nucleation  Process 


The  relative  intensities  of  the  features  of  various  cluster 
geometries  can  reveal  something  of  the  nucleation  processes  in  the 
molecular  jet.  The  Ben(S)2  data  are  particularly  useful  in  focusing 
on  these  processes.  In  general,  two  geometries  can  exist  for  this 
complex:  a  solvent  molecule  can  be  on  either  side  of  the  ring 
(isotropic  cluster)  or  both  solvent  molecules  can  be  on  the  same  side 
of  the  ring  (anisotropic  cluster).  Although  homogeneous  nucleation 
can  conceivably  create  either  type  of  cluster,  inhomogeneous 
nucleation  can  only  create  anisotropic  clusters  of  this  size.  By 
examining  the  relative  intensities  of  the  spectral  features  of  each 
type  of  cluster,  one  can  determine  the  relative  concentration  of 
isotropic  vs.  anisotropic  clusters  in  the  jet.  In  this  way  the 
nucleation  processes  in  the  jet  can  be  elucidated. 

The  spectra  of  Ben(C2Hg)2  «nd  Ben(C^Hg)2  (figures  4  and 
8)  show  that  there  are  virtually  no  isotropic  clusters  of  these 
species.  Since  homogeneous  nucleation  should  produce  both  isotropic 
and  anisotropic  clusters,  inhomogeneous  nucleation  appears  to  be  the 
dominant  nucleation  process  for  these  systems  under  the  present  beam 
conditions.  By  contrast,  the  Ben(CH^)2  spectrum  (figure  l)  shows 
both  isotropic  and  anisotropic  clusters  in  roughly  equal  numbers. 
Therefore,  homogeneous  nucleation  does  occur  for  the  benzene-methane 
system. 

The  lack  of  homogeneous  nucleation  in  Ben(C2Hg)2  and 
Ben(C_Ho)o  indicates  that  the  number  of  solvent  monomers  relative 

j  o 

to  dimers  and  larger  solvent  clusters  is  probably  small  for  ethane  and 
propane  solvents.  However,  based  on  the  calculated  binding  energies 
the  relative  number  of  methane  monomers  should  also  be  small. 
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Therefore,  it  is  unlikely  that  enough  methane  monomers  exist  in  the 
beam  to  explain  the  amount  of  homogeneous  nucleation  observed  for 
Ben(CH^)2. 

It  is  difficult  to  get  data  on  the  relative  concentrations  of 

molecular  clusters  in  a  jet  expansion,  especially  for  the  hydrocarbon 

systems.  Some  data  are  available  on  the  concentration  of  dimers 

relative  to  monomers  for  ethylene  under  rather  modest  expansion 

conditions.  (Pq=  5-10  atm,  d  =  35  ym  orifice).  The  apparent 

ratio  of  monomers  to  dimers  is  roughly  4:1  for  these  conditions:  The 

2 

relative  concentration  of  clusters  scales  with  (P  d)  and  thus  our 

o 

expansion  conditions  will  produce  a  factor  of  15  higher  concentration 
of  clusters.  The  observed  4:1  ratio  for  ethylene  is,  however,  heavily 
distorted  by  fragmentation  caused  by  the  electron  impact  ionization 
technique:  for  example,  almost  all  observed  dimers  arise  from  higher 
clusters.  Hence,  the  actual  ratio  of  monomer  to  dimer  neutrals  is  not 
well  estimated  by  those  methods;  dimer  and  higher  clusters  could 
easily  dominate  the  distribution  of  species  in  their  expansion  and 
most  certainly  dominate  under  our  more  extreme  expansion  conditions. 

If  the  homogeneous  nucleation  of  BenCCH^)^  is  not  occuring 
through  the  interactions  of  methane  monomers  and  benzene  (and  a  third 
body  to  carry  off  the  binding  energy)^ this  nucleation  must  be 
occuring  through  the  interaction  of  dimers  (or  larger  clusters)  with 
benzene.  Upon  interaction  with  benzene,  the  methane  dimer  can 
dissipate  part  of  the  binding  energy  through  vibrational 
predissociation  (VP).  This  interactive  collision  leaves  one  methane 
bound  to  benzene  while  the  other  methane  carries  off  a  portion  of  the 
excess  binding  energy. 

The  binding  energies  of  the  solvent  dimers,  as  well  as  those 
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for  benzene  with  the  three  solvents,  are  given  in  Table  XI.  The 
values  are  derived  from  computer  modeling  of  the  clusters.  The  data 
suggest  that  only  about  half  of  the  binding  energy  of  a  solvent  to 
benzene  can  be  dissipated  through  VP  of  the  solvent  dimer. 

If  the  binding  energy  does  not  cause  VP  of  thu  solvent  dimer 
before  a  third  body  interacts  with  the  complex  and  carries  away  the 
excess  binding  energy,  the  resultant  clusters  will  be  anisotropic. 
Therefore,  the  nucleation  of  benzene  with  ethane  and  propane  is 
inhomogeneous  because  the  solvent  dimers  do  not  undergo  VP  quickly 
enough  during  the  collision/ formation  process. 

Based  on  the  above  reasoning,  one  would  expect  that  a  large 
solute-solvent  binding  energy  would  result  in  an  increased  number  of 
homogenously  nucleated  (.isotropic)  clusters.  Evidence  is  presented  in 
II  for  the  toluene-solvent  systems  which  support  this  conclusion. 

V.  CONCLUSIONS 

Two-color  TOFMS  studies  have  been  employed  to  determine  the 
general  geometry  and  symmetry  of  small  benzene-hydrocarbon  clusters: 
comparison  of  cluster  0^  and  transitions  makes  possible  the 
experimental  identification  of  cluster  geometry  (symmetry)  with 
spectral  shift  and  vibrational  features.  Through  computer  modeling  of 
van  der  Waals  clusters  of  benzene  and  small  hydrocarbons,  a 
correlation  between  detailed  cluster  geometry  and  spectral  features 
has  proven  possible.  The  specific  assignments  lead  to  a  general 
pattern  for  the  relationship  between  geometry  and  spectral  shift.  The 
spectral  shift  in  these  systems  seems  to  be  largely  due  to  the 
interaction  of  the  solvent  with  the  aromatic  w -electron  system  of 


benzene.  This  would  explain  the  observations  that  spectral  shift  and 
calculated  binding  energy  are  not  directly  related.  Of  course,  the 
details  of  these  calculations  are  a  function  of  the  model  chosen  for 
the  potential  (atom  -  atom  exp-6). 

The  assignment  of  geometry  to  clusters  with  a  specific  spectral 
shift  has  also  made  possible  a  better  understanding  of  the  nucleation 
processes  in  the  jet.  We  suggest  that  due  to  the  relatively  low 
concentration  of  solvent  monomers  in  the  system,  homogeneous 
nucleation  can  occur  through  the  VP  of  solvent  dimers  as  the  dimers 
interact  with  benzene.  If  the  binding  energy  of  the  solvent  dimer  is 
not  small  relative  to  the  solute-solvent  binding  energy,  VP  will  not 
occur  and  anisotropic  geometries  of  the  solute-solvent  clusters  will 
comprise  the  major  cluster  species  in  the  beam. 
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Aliphatic -Aliphatic^ 


c-c 

131096 

8890353 

3.421 

C-H 

47830 

7562708 

3.94 

H-H 

15028 

6390528 

4.643 

Aliphatic-Aromatic^ 


C-C 

162645 

16520286 

3.5105 

C-H 

47864 

6434085 

3.805 

H-H 

13071 

2497533 

4.1915 

Aromatic-Aromatic*^ 


C-C 

201786 

30698427 

3.6 

C-H 

47897 

5473891 

3.67 

H-H 

11368 

976084 

3.74 

a) 

b) 


From  the  methane-methane  parameters  of  ref.  11 

The  parameters  are  found  by  combining  the  benzene-benzene  and  methane-methane 
parameters  as  follows: 


Sm  ^^Sb 


From  the  benzene-benzene  parameters  of  ref.  10 


TABLE  II 


The  prominent  peaks  in  the  spectra  of  Ben(CH^)j^  and  Ben(CH^)2  in  the  region  of 
6q  (see  figure  2). 


Species 

Energy 
(vac  cm 

Energy  relative 

To  Ben  6^  (cm“^) 

Energy  relative 

To  cliister  6g 

Assignment^ 

B€n(CH^)^ 

38567.6 

-41 

0 

38594.6 

27.3 

^0^0 

38599.9 

32.3 

^yo 

38619.0 

51.4 

BenCCH^)^ 

38527.1 

-81.4 

0 

I  6^ 

38557.1 

30.6 

38568.1 

-40.4 

0 

M 

M 

o^ 

O  (-• 

TABLE  III 


The  atcroic  positions  for  the  minimum  energy  configurations  of 
Ben(CH^)j^  and  Ben(CH^)2  (see  figure  3).  Coordinates  are  in  &  with  the  origin 
at  the  center  of  benzene,  the  X  axis  passes  through  a  carbon  atom  on  the 
benzene,  the  Y  axis  bisects  a  C-C  bond  on  benzene  and  the  Z  axis  is  coincident 
with  the  6-fold  synmetry  axis  of  benzene. 


Atom 


X(&) 


y(8) 


Z(&) 


Energy 

(cm"^) 


TABLE  IV 

The  prOToinent  peaks  in  the  spectra  of  Ben(C2Hg)2^  and  Ben(C2Hg)2  in  the  region  of 
6q  (see  figure  4) . 


Species 

Energy 

Energy  relative 

Energy  relative 

Assignment 

(vac  cm~^) 

to  Ben  6Q(cm 

to  cluster  6j(cm 

TABLE  V 


The  prominent  peaks  in  the  spectra  of  Ben(C_H^)<  and 

<1  o  1 

Ben(C2Hg)2  in  the  0^  region  (see  figure  5). 


Species 

Frequency 

Energy  relative 

Energy  relative 

Assignment® 

( vac  cm~^ ) 

to  Ben(0Q(cm 

to  cluster  0^ 

TABLE  VI 

The  atomic  positions  for  the  nrimnum  energy  configurations  of  Ben(C2H^)]^ 
(see  figure  6).  The  benzene  coordinates  are  the  same  as  in  Table  III. 


Atom  X(X)  Y(R)  Z(X)  Energy 

(cm"^) 


TABLE  VII 


The  atxxnic  positions  for  the  ndniiaim  energy  configurations  of  Ben(C2H^)2 
(see  figure  7) .  The  benzene  coordinates  are  the  same  as  in  Table  III . 


Atom  X(X)  Y(8)  Z(S)  Energy 

(cm'^) 


Configuration  I 

Ligand  1  Cl 

0 

C2 

0 

ffl. 

-  .8989 

H2 

.8989 

H3 

0 

H4 

-  .8989 

H5 

0 

H6 

.8989 

Ligand  2  Cl 

.001 

C2 

.001 

HI 

-  .8979 

H2 

.8999 

H3 

.001 

H4 

.001 

H5 

-  .8999 

H6 

.8979 

Configuration  II 

Ligand  1 

Cl 

-  2.624 

C2 

-  3.177 

HI 

-  3.318 

H2 

-  1.643 

H3 

-  2.516 

H4 

-  4.257 

H5 

-  3.284 

H6 

-  2.482 

Ligand  2 

Cl 

.374 

C2 

.771 

HI 

.170 

H2 

1.196 

H3 

-  .528 

H4 

-  .051 

H5 

1.674 

H6 

.976 

-988 


0 

3.396 

0 

4.926 

-  .519 

3.031 

-  .519 

3.031 

1.037 

3.031 

.519 

5.291 

-  1.037 

5.291 

.519 

5.291 

0 

8.346 

0 

9.876 

-  .518 

7.981 

-  .518 

7.981 

1.038 

7.981 

-  1.038 

10.241 

.518 

10.241 

.518 

10.241 

-1512 


-  1.486 

3.000 

-  1.978 

4.339 

-  1.766 

2.193 

-  1.948 

2.818 

-  .392 

3.029 

-  1.514 

4.527 

-  3.072 

4.310 

-  1.698 

5.145 

-  .391 

4.808 

.268 

3.485 

-  1.458 

4.636 

-  .286 

5.530 

.098 

5.203 

.163 

2.763 

-  .220 

3.091 

1.336 

3.657 
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TABLE  VIII 


The  prominent  peaks  in  the  spectnin  of  Ben(C^)^  and  Ben(C2Hg)2  in  the  6q 
region  (see  figure  8). 


TABLE  X 


The  prorament  peaks  in  the  spectum  of  Ben(C2Hg)2  in  the  6q  and  Oq  regions 
(see  figure  10) . 


Energy 
(vac  cm~^) 


Energy  relative 
to  Ben  feature 
(cm"b 


Energy  relative 
to  cluster  feature 

(cm"b 


Assignment 


TABLE  XI 

The  ratio  of  the  solvent  dimer  binding  energy  to  the  benzene -solvent  binding  energy  (BE) 


Solvent 

Solvent  dimer  B.E.^ 

Benzene-  so  Ivent 

Ra 

(cm~^) 

B.E.  (cm 

methane 

ethane 

propane 


277 

599 

.47 

503 

778 

.65 

789 

1044 

.75 

All  values  are  calculated  by  computer  modeling  of  the  complex 
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FIGURE  CAPTIOUS 


FIGURE  1 

Mass  selective  absorption  spectra  of  Ben(CH^)^,  BenlCH^J^ 
BeniCH^J^  in  the  region  of  the  6^  transition  of  benzene.  The 
energy  scale  is  relative  to  the  Ben  6^  transition  (58608.5  cm~^). 
PeaK  positions  are  given  in  Table  II.  The  peak  positions  for 
Ben(CH^)^  are  not  tabulated. 

FIGURE  2 

Mass  selective  absorption  spectrum  of  BenCCH^)^  in  region  of 
Oq  transition  of  benzene.  The  energy  scale  is  relative  to  the 
Ben(0Q  (58086.1  cm  ^ ^ .  The  presence  of  an  origin  for 
configuration  II  of  Ben(CH^)2  (ca.  >8042  cm~^)  indicates  that  the 
5-fold  axis  of  benzene  is  destroyed  in  this  geomeiry. 


FIGURE  5 

Minimum  energy  configurations  of  Ben;.CH^^^  and  Ben(,CH^)2.  Atomic 
positions  are  given  in  Table  ill.  Note  that  the  5-fold  axis  of 
benzene  is  preserved  in  Ben(CH^)^  and  geometry  I  of  Ben(CH^)2.  bu 
not  in  geometry  II  of  Ben(CH^)^. 


FIGURE  4 


Maas  selective  absorption  spectra  of  BenCC^Hg)^  and  Ben(C2Hg)2 
in  the  6q  region  of  benzene.  The  energy  scale  is  relative  to  the 
Ben  6q  (58608.5  cm”^ ) .  The  peak  position  indicated  are  tabulated 
in  Table  IV. 

FIGURE  5 

Mass  selective  absorption  spectra  of  BenCC^Hg)^  and  Ben(C2Hg)2 
in  the  Oq  region  of  benzene.  The  energy  scale  is  relative  to  the 
Ben  Oq  (58086.1  cm  ^ ) .  Note  that  only  configuration  II  exhibits 
an  origin  and  vdW  bends.  Peak  positions  are  given  in  Table  V. 

FIGURE  6 

Minimum  energy  geometries  of  Ben(C2Hg)^.  Atomic  positions  are 
given  in  Table  VI.  Note  that  in  geometry  I  the  5-fold  axis  of  benzene 
is  preserved,  but  in  II  it  is  not. 

FIGURE  7 

Minimum  energy  geometries  of  Ben(C2Kg)2.  Atomic  positions  are 
given  in  Table  VII.  Note  that  geometry  I  preserves  the  5-fold  symmetry 
of  benzene,  but  geometry  II  dof’s  not. 

FIGURE  8 

Mass  selective  absorption  spectra  of  Ben(C^Hg)^  and  Ben(C^Hg)2 
in  the  6^  region  of  benzene.  The  energy  scale  is  relative  to  the 
Ben  6g(58608.5  cm~M.  The  peak  positions  indicated  are  tabulated 
in  Table  VIII. 


FIGURE  y 


Minimum  energy  geometries  of  Ben(.C^Hy;^  and  Ben(,C^Hy)^.  Atomic 
positions  are  given  in  Table  IX.  Note  that  the  3-fold  symmetry  of 
benzene  is  nearly  preserved  by  the  atoms  closest  to  benzene  for  both 
species. 

FIGURE  10 

Mass  selective  absorption  spectrum  of  Ben\.C^Hyj.^  in  the  and 
Oy  regions.  The  energy  scales  are  relative  to  the  Ben  Oq  (3tibU6.3 
cm~^ )  and  the  Ben  0^  (jSUtib.l  cm”^;.  Peak  positions  aru  given 
in  Table  X.  Note  that  both  geometries  of  BeniC.^H^)^  destroy  the 
3-fold  symmetry  of  benzene. 
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